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The set of Bell-diagonal states for two qubits can be depicted as a tetrahedron in three dimensions. 
We consider the level surfaces of entanglement and quantum discord for Bell-diagonal states. This 
provides a complete picture of the structure of entanglement and discord for this simple case and, in 
particular, of their nonanalytic behavior under decoherence. The pictorial approach also indicates 
how to show that discord is neither convex nor concave. 

PACS numbers: 03.67.-a, 03.67.Mn, 03.65.Ud, 03.65.Yz 
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Maintenance of quantum coherence is clearly im- 
portant for quantum-information-processing protocols. 
Noise and decoherence, by turning pure states into mixed 
states, generally destroy quantum coherence. Efficient 
representation of quantum information requires that a 
quantum-information-processing system be composed of 
parts For multi-partite systems, quantum coherence 
is related to nonclassical correlations between the parts. 

One kind of nonclassical correlation is entangle- 
ment [2j. A pure quantum state is unentangled if it is 
a product of pure states for each part. A mixed state is 
unentangled (separable) if it can be written as an ensem- 
ble of such product states. Entanglement is the crucial 
resource for such quantum-information-processing pro- 
tocols as quantum key distribution, teleportation, and 
super-dense coding 0. 

Operational measures of entanglement are notoriously 
difficult to calculate for mixed states; even the bound- 
ary between separability and entanglement is difficult to 
characterize. One can say, however, that the set of sepa- 
rable states is a convex set, is invariant under local uni- 
tary operations, and has dimension as large as the space 
of mixed states f2| . 

Separable states have nonzero measure in the space of 
all states 0] . In a decoherence process that involves decay 
to a separable equilibrium state that does not lie on the 
boundary between separability and entanglement, the de- 
cohering state will cross that boundary before reaching 
the equilibrium state. This phenomenon, dubbed "sud- 
den death of entanglement" 0, [f| , is the generic expec- 
tation in view of the geometry of separable states. 

Separable states can have nonclassical correlations 
even though they are unentangled. A state with only 
classical correlations, often called a classical state, is one 
that is diagonal in a product basis, for then the corre- 
lations are described by a joint probability distribution 
for classical variables of the parts. These purely classical 
states are a set of measure zero, as is suggested by the fact 
that any classical state can be perturbed infinitesimally 
to become nonclassical by making two of the eigenvectors 
infinitesimally entangled and is proved rigorously in Q . 

A variety of measures have been proposed to quantify 



nonclassical correlations for bipartite systems m 
ways that can be nonzero for separable, but nonclassi- 
cal states. Nonclassical, but perhaps separable correla- 
tions have been related to exponential speed-ups in the 
"power-of-one-qubit" model [1 01 ] of mixed-state quantum 
computation [111 ], but the relation remains tenuous (12J. 

One can use decoherence mechanisms to explore the 
nooks and crannies of nonclassical-correlation measures. 
There is no sudden death [fj, as is suggested by the ab- 
sence of open sets of classical states, but the nonanalyt- 
icity of nonclassical measures points to the possibility of 
sudden changes in derivatives. Investigation of the be- 
havior of nonclassical measures under decoherence has 
begun 

SEMI, 

with a focus on the action of decoher- 
ence within the class of two-qubit states that are diagonal 
in the Bell basis. This focus is motivated by the fact that 
entanglement measures and nonclassical-correlation mea- 
sures can be calculated explicitly for the Bell-diagonal 
states, thus allowing one to determine how these mea- 
sures change under decoherence. 

The Bell-diagonal states are a three-parameter set, 
whose geometry, including the subsets of separable and 
classical subsets, can be depicted in three dimensions [H, 
[lfl ] . Level surfaces of entanglement and nonclassical mea- 
sures can be plotted directly on this three-dimensional 
geometry. The result is a complete picture, for this sim- 
ple case, of the structure of entanglement and nonclas- 
sicality. We suggest that it is more illuminating to use 
this picture to explain how measures of entanglement and 
nonclassicality change along the one-dimensional trajec- 
tories traced out by decohering states, rather than the 
other way around. Hence we review and expand the pic- 
torial approach here. 

The Bell-diagonal states of two qubits, A and B, have 
density operators of the form 

PAB = \(l + J2cj af ® trf\=^\ab\Pab){Pab\ , (1) 
^ 3=1 ' a,b 

where the aj's are Pauli operators. The eigenstates are 
the four Bell states \/3 ab ) = (\0,b) + (-1)°|1, 1 © 6))/V2, 
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FIG. 1: Geometry of Bell-diagonal states. The tetrahedron 
T is the set of valid Bell-diagonal states. The Bell states 
\Pab) sit at the four vertices, the extreme points of T. The 
green octahedron O, specified by |ci| + |c2 1 + |C3 1 < 1 (\ a b < 
1/2), is the set of separable Bell-diagonal states. There are 
four entangled regions outside O, one for each vertex of T, in 
each of which the biggest eigenvalue X a b is the one associated 
with the Bell state at the vertex. Classical states, i.e., those 
diagonal in a product basis, lie on the Cartesian axes. 



with eigenvalues 



Xab = \(}+ (-l) Q Cl - (-l) a+b C 2 + (-1) 6 C 3 ) 



(2) 



Any two-qubit state satisfying (af) = = (o~f), 
i.e., having maximally mixed marginal density operators 
PA = 1/2 = ps, can be brought to Bell-diagonal form 



by using local unitary operations on the two qubits to 
diagonalize the correlation matrix (a^ (g) crj?)- 

A Bell-diagonal state is specified by a 3-tuplc 
(ci,C2,C3). The density operator pab must be a posi- 
tive operator, i.e., X a b > 0; the resulting region of Bell- 
diagonal states is the state tetrahedron T in Fig. [T] Sep- 
arable Bell-diagonal states are those with positive partial 
transpose |2j. Partial transposition changes the sign of 
C2, so operators with positive partial transpose occupy 
the reflection of T through the plane c 2 = 0; the region 
of separable Bell-diagonal states is the intersection of the 
two tetrahedra, which is the octahedron O of Fig. [Tl [lij]. 

The entanglement of formation £ 0, 17| is a mono- 
tonically increasing function of Wootters's concur- 
rence C [13, which for Bell-diagonal states, is given by 
C = max(0, 2A max — 1), where A max = max A f,. The con- 
currence and the entanglement of formation are convex 
functions on T. They are zero for the separable states in 
the octahedron O. In each of the four entangled regions 
outside O, C and £ are constant on planes parallel to the 
bounding face of O and increase as one moves outward 
through these planes toward the Bell-state vertex. 

Quantum discord was introduced by Ollivier and 



Zurek Q ■ We restrict attention to it because of its promi- 
nence among measures of nonclassical correlations and 
because it has been a focus of recent work on decoher- 
ence and nonclassical correlations SEMI. 

To define quantum discord, one starts with the quan- 
tum mutual information, I = S(pa) + S{pb) — S(pab) — 
S(pb) — S(B\A), where S(p) = — tr(plog 2 p) is the von 
Neumann entropy of p and S(B\A) = S(pab) — S{pa) 
is a conditional quantum entropy. The quantum mutual 
information is regarded as quantifying the total correla- 
tions in the joint state pab- 

The quantum mutual information of Bell-diagonal 
states, 



1=2- S(pab) = Xab l0 S2(^ab) 



(3) 



a . b 



is a convex function on T. It has smooth level surfaces 
that bulge outward toward the vertices of T. 

The next step is to quantify purely classical correla- 
tions in terms of information from measurements. One 
imagines measuring on A a POVM consisting of rank- 
one POVM elements E k = Dq k \k)(k\ [Hj], where D is 
the dimension of A and the q k make up a normalized 
probability distribution. The probability to get result k 
is pk = Dqk(k\pA\k), and the post-measurement state of 
B is ps\k = (k\PAB\k) / (k\pA\k) . Minimizing the average 
entropy of B, given result k, over measurements on A, 
yields a classical conditional entropy 



S(B\A) = mmy2p k S(p B \k) I 



(4) 



minimizing chooses the measurement of A that extracts 
as much information as possible about B. The corre- 
sponding mutual-information-like quantity C = S(ps) — 
S(B\A) is the measure of classical correlations. 
For Bell-diagonal states, we have 



C = l-H 2 



1 



1 



c, , , 1 
-log 2 (l + c) + - 



log 2 (l-c) , 



(5) 

(1 — p) log 9 (l — p) is the bi- 
This C. a convex 
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where H2 (p) = —p log 2 p - 
nary entropy and c = max 
function on T, is constant on the surfaces of cubes (or 
the portion of such a cube in T) centered at the origin of 
Fig. [TJ — this introduces nonanalyticity — and C increases 
monotonically with the size of the cube. 

Discord is defined as the difference of I and C, 



V =1 - C = S(B\A) - S(B\A) 



(6) 



thus capturing a notion of nonclassical correlations. 
Since C is generally asymmetric between A and B, so 
also is the discord; this means, in particular, that dis- 
cord, as defined, vanishes if and only if pab is diagonal 
in a conditional product basis \ef) ® |/^.), rather than 
only in a product basis \ef) (g> \ f k )- Bell-diagonal states 
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being symmetric between A and B, however, discord is 
zero only for classical states, which lie on the Cartesian 
axes in Fig. [TJ 0- 

Figure [5]plots level surfaces of discord for Bell-diagonal 
states. From these plots, it is clear that discord is quite a 
different beast from entanglement of formation, quantum 
mutual information, and the measure of classical corre- 
lations. Whereas £, I, and C generally increase outward 
from the origin, T> increases away from the Cartesian 
axes, capturin g an entropic notion of distance from clas- 
sical states Il2|. In particular, as one moves outward 
along one of the constant-discord tubes of Fig. [2j the clas- 
sical correlations and the total correlations of the quan- 
tum mutual information increase, but their difference, 
the nonclassical correlations as measured by discord, re- 
mains constant. At the vertices of O, X = C = 1 and 
V = £ = 0. At the Bell-state vertices of T, T = 2 
and C = T> = £ = 1, this being the maximum value 
of discord for two qubits. In addition, £ , I, and C are 
all convex, whereas discord is neither concave nor con- 
vex, as is evident from the plots in Fig. [5] one can mix 
two positive-discord states to get a zero-discord classical 
state, and one can mix two zero-discord classical states 
on different axes to get a positive-discord state [ilf . 

Mazzola, Piilo, and Maniscalco Iff recently investi- 
gated the dynamics of classical and nonclassical correla- 
tions, as measured by discord, for two qubits under deco- 
herence processes that preserve Bell-diagonal states. In 
particular, they considered independent phase-flip chan- 
nels for the two qubits. The phase flips are implemented 
mathematically by random applications of a z operators 
to the qubits. This decoherence process leaves C3 un- 
changed, but flips the signs of C\ and c-i randomly, lead- 
ing to exponential decay of c\ and ci at the same rate. 
Mazzola and collaborators found that for the initial con- 
ditions they considered, the entanglement of formation 
decays to zero in a finite time — sudden death of entan- 
glement [IJ — but that the discord remains constant for 
a finite time and then decays, reaching zero at infinite 
time. This situation is depicted in terms of the surfaces 
of constant discord in Figure [3J The decohering-state 
trajectory is a straight line that runs along a tube of 
constant discord, until it encounters an intersecting tube, 
after which the discord decreases to zero when the state 
becomes fully classical. 

This behavior is generic for flip channels and initial 
conditions on edges of the state tetrahedron. We focus 
here on the phase-flip channel with initial conditions in 
the (+, — , +)-octant, but analogous considerations apply 
to the other flip channels (bit and bit-phase) and to initial 
conditions on the other edges of T. Consider then initial 
conditions anywhere along the edge of T in this octant: 
ci(0) = 1 and < -c 2 (0) = c 3 (0) < 1. The trajectory 
under phase flips is a straight line C3 = 03(0) = — C2/C1. 
Along this straight line, the eigenvalues \ a b factor into 
products of probabilities, (1 ± ci)/2 and (1 ± C3)/2, thus 
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FIG. 2: Surfaces of constant discord: (a) V = 0.03; (b) V = 
0.15; (c) T> = 0.35. The level surfaces consist of three inter- 
secting "tubes" running along the three Cartesian axes. The 
tubes are cut off by the state tetrahedron T at their ends, 
and they are squeezed and twisted so that at their ends, they 
align with an edge of T. As discord decreases, the tubes col- 
lapse to the Cartesian axes [l2l ]. As discord increases, the 
tube structure is obscured, as in (c): the main body of each 
tube is cut off by T; all that remains are the tips, which reach 
out toward the Bell-state vertices. 
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FIG. 3: Trajectory (red) of a Bell-diagonal state under ran- 
dom phase flips of the two qubits; initial conditions are 
ci(0) = 1, -c 2 (0) = c 3 (0) = 0.3. The trajectory is the 
straight line c$ = cs(0) = 0.3 = —Qz/ci. For clarity, only 
the (+, — , +)-octant is shown. A constant-discord surface is 
plotted for the discord value of the initial state. Faces of the 
yellow state tetrahedron T and the green separable octahe- 
dron O are also shown. The straight-line trajectory proceeds 
along a tube of constant discord till it encounters the vertical 
tube at ci = 0.3, after which discord decreases monotonically 
to zero when the trajectory reaches the C3 axis. Entangle- 
ment of formation decreases monotonically to zero when the 
trajectory enters O at ci = 0.7/1.3 = 0.54. 



making S(pab) the entropy of two independent binary 
random variables with these probabilities. This yields a 
quantum mutual information X = 2 — i?2[(l + c 3)/2] — 
H2KI + ci)/2]. Furthermore, along the trajectory c = 
max(ci,C3). The result is that the trajectory initially 
runs along a tube of constant discord 



V = 1 - Hi 



l + c 3 



(7) 



for ci > C3. When c\ — C3, the trajectory encounters 
another tube, after which, for c\ < C3, the discord de- 
creases monotonically as V = 1 — H 2 [(l + c 1 )/2] as c\ 
decreases. Meanwhile, the entanglement of formation de- 
creases monotonically from its initial value to a sudden 
death at c\ = (1 — c 3 )/(l + c 3 ) 



The situation investigated in [151 ] is surely interesting: 
under decoherence, nonclassical correlations remain con- 
stant for a finite time interval. This situation is, how- 
ever, a special one, as can be seen from the surfaces 
of constant discord; the trajectories considered here are 
the only straight lines in parameter space that stay on 
a surface of constant discord. Indeed, the pictorial ap- 
proach can provide a complete understanding of how en- 
tanglement and nonclassicality change under decoherence 
within the set of Bell-diagonal states. 
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